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-Abstract- 

In this paper we study similarity measures for moving curves which can, for example, model 
changing coasthnes or retreating glacier termini. Points on a moving curve have two parameters, 
namely the position along the curve as well as time. We therefore focus on similarity measures for 
surfaces, specihcally the Frechet distance between surfaces. While the Frechet distance between 
surfaces is not even known to be computable, we show for variants arising in the context of moving 
curves that they are polynomial-time solvable or NP-complete depending on the restrictions 
imposed on how the moving curves are matched. We achieve the polynomial-time solutions by 
a novel approach for computing a surface in the so-called free-space diagram based on max-flow 
min-cut duality. 
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Introduction 


Over the past years the availability of devices that can be used to track moving objects has 
increased dramatically, leading to an explosive growth in movement data. Naturally the 
goal is not only to track objects but also to extract information from the resulting data. 
Consequently recent years have seen a significant increase in the development of methods 
extracting knowledge from moving object data. 

Tracking an object gives rise to data describing its movement. Often the scale at which 
the tracking takes place is such that the tracked objects can be viewed as point objects. Cars 
driving on a highway, birds foraging for food, or humans walking in a pedestrian zone: for 
many analysis tasks it is sufficient to consider objects as moving points. Hence the most 
common data sets used in movement data processing are so-called trajectories: sequences of 
time-stamped points. 

However, not all moving objects can be reasonably represented as points. A hurricane can 
be represented by the position of its eye, but a more accurate description is as a 2-dimensional 
region which represents the hurricane’s extent. When studying shifting coastlines, reducing 
the coastline to a point is obviously unwanted: one is actually interested in how the whole 
coast line moves and changes shape over time. The same holds true when studying the 
terminus of a glacier. In such cases, the moving object is best represented as a polyline rather 
than by a single point. In this paper we hence go beyond the basic setting of moving point 
objects and study moving complex, non-point objects. Specifically, we focus on similarity 
measures for moving curves, based on the Frechet distance. 


Definitions and Notation. The Frechet distance is a well-studied distance measure for shapes, 
and is commonly used to determine the similarity between two curves A and B : [0,1] —>■ M". 
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Computing the Similarity Between Moving Curves 


A natural generalization to more complex shapes uses the definition of Equation where the 
shapes A and B have type X —> M". 

Dm(A,B)= inf sup||A(x)-B(/r(a:))|| (1) 

Here, || • || : K" —>■ M is a norm such as the Euclidean 
norm (L^) or the Manhattan norm {L^). The match¬ 
ing ^ ranges over orientation-preserving homeomorph- 
isms (possibly with additional constraints) between 
the parameter spaces of the shapes compared; as such, 
it defines a correspondence between the points of the 
compared shapes. A matching between surfaces with 
parameters p and t is illustrated in Figure Given 
one such matching we obtain a distance between A 
and B by taking the largest distance between any two 
corresponding points of A and B. The Frechet distance 
is the infimum of these distances taken over all possible matchings. For moving points or 
static curves, we have as parameter space X = [0,1] and for moving curves or static surfaces, 
we have X = [0,1]^. We can define various similarity measures between shapes by imposing 
further restrictions on p,. 

In practice a curve is generally represented by a sequence of P -I- 1 points. Assuming a 
linear interpolation between consecutive points, this results in a polyline with P segments. 
Analogously, a moving curve is a sequence of T -|- 1 polylines, each of P segments. We also 
interpolate the polylines linearly, yielding a bilinear interpolation, or a quadrilateral mesh 
of P X T quadrilaterals. 

Related Work. The Frechet distance or related measures are frequently used to evaluate the 
similarity between point trajectories [SI Ellin]. The Frechet distance is also used to match 
point trajectories to a street network mm- The Frechet distance between polygonal curves 
can be computed in near-quadratic time 13 El 13 EH], and approximation algorithms mm 
have been studied. 

The natural generalization to moving (parameterized) curves is to interpret the curves 
as surfaces parameterized over time and over the curve parameter. The Frechet distance 
between surfaces is NP-hard m, even for terrains m- In terms of positive algorithmic 
results for general surfaces the Frechet distance is only known to be semi-computable mm- 
Polynomial-time algorithms have been given for the so called weak Frechet distance [T] and 
for the Frechet distance between simple polygons m and so called folded polygons m- 

When interpreting moving curves as surfaces it is important to take the different roles of 
the two parameters into account: the first is inherently linked to time and the other to space. 
This naturally leads to restricted versions of the Frechet distance of surfaces. For curves, 
restricted versions of the Frechet distance were considered ElES- For surfaces we are not 
aware of similar results. 

1.1 Results 

We refine the Frechet distance between surfaces to meaningfully compare moving curves. To 
do so, we restrict matchings to be one of several suitable classes. Representative matchings 
for the considered classes together with the running times of our results are illustrated in 
Figure]^ 


A Bop 




M Figure 1 A matching p between sur¬ 
faces A and B drawn as a homeomorph- 
ism between their parameter spaces. 
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The simplest class of matchings consists of a single pre¬ 
defined identity matching yi{p,t) = {p,t). Hence, to compute 
the identity Frechet distance, we need only determine a pair of 
matched points that are furthest apart. It turns out that one of 
the points of a furthest pair is a vertex of a moving curve (i.e. 
quadrilateral mesh), allowing computation in 0{PT) time, see 
Section [21 

We discuss the synchronous constant Frechet distance in Sec¬ 
tion Here we assume that the matching of timestamps is known 
in advance, and the matching of positions is the same for each 
timestamp, so it remains constant. Our algorithm computes the 
positional matching minimizing the Frechet distance. 

The synchronous dynamic Frechet distance considered in Sec- 
tionj^also assumes a predefined matching of timestamps, but does 
not have the constraint of the synchronous constant class that 
the matching of positions remains constant over time. Instead, 
the positional matching may change continuously over time. 

Finally, in Section we consider several cases where neither 
positional nor temporal matchings are predefined. The three 
considered cases are the asynchronous constant, asynchronous 
dynamic, and orientation-preserving Frechet distance. The asyn¬ 
chronous constant class of matchings consists of a constant (but 
not predefined) matching of positions, as well as timestamps 
whereas in the asynchronous dynamic class of matchings, the 
positional matching may change continuously. In the orientation¬ 
preserving class, matchings range over orientation preserving 
homeomorphisms between parameter spaces, given that the 
corners of the parameter spaces are aligned. 

The last three classes are quite complex, and we give construc¬ 
tions proving that approximating the Frechet distance within a 
factor 1.5 is NP-hard under these classes. For the asynchronous 
constant and asynchronous dynamic classes of matchings, this 
result holds even for moving curves embedded in whereas the 
result for the orientation-preserving case holds for embeddings 
in 


Although we do not discuss classes where positional match¬ 
ings are known in advance, these symmetric variants can be 
obtained by interchanging the time and position parameters for 
the discussed classes. Deciding which variant is appropriate for 
comparing two moving curves depends largely on how the data is 
obtained, as well as the use case for the comparison. For instance, 
the synchronous constant variant may be used on a sequence 
of satellite images which have associated timestamps. The syn¬ 
chronous dynamic Frechet distance is better suited for sensors 
with different sampling frequencies, placed on curve-like moving 
objects. 



Identity 

0{PT) 



Synchronous Constant 
0(P^riog(Pr)) 



Synchronous Dynamic 
0(P^riogPlog(PT)) 



Asynchronous Constant 
NP-complete 



Asynchronous Dynamic 
NP-hard 



Orientation-Preserving 

NP-hard 


H Figure 2 The time com¬ 
plexities of the considered 
classes of matchings. 
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\~ 2 ] Identity Matchings 


Suppose we are given a single predefined matching ^ between the moving 
curves A and B : [0, P] x [0, T] —> K". We can compute the Frechet distance 
under this matching if we can find the corresponding points of A and B that 
are furthest apart, see Equation 

D^{A,B)= sup ||A(x) - P(Ai(a;))|| (2) 

a;G[0,P]x[0,T] 


If A and B are both of size P xT, then the identity matching /i(p, t) = {p, t) allows us to 
simplify A{x) — B(p{x)) into C{x) = A{x) — B{x) where C is again a quadrilateral mesh 
of size P X T. The Frechet distance for the identity matching depends only on the point 
on C that is furthest from the origin. To see this, consider a single quadrilateral, the point 
furthest from the origin must be one its four corners since all points on the quadrilateral lie 
within the convex hull of its four corner points. Hence, it suffices to check only the distance 
to the origin for the 0{PT) vertices of C, see Equation]^ The Frechet distance under the 
identity matching can then be computed in 0{PT) time. 


Did{A,B)= sup II H(a:) - P(a;) II = sup ||C'(a;)|| (3) 

a:G[0.P]x[0,T] '-v-^ a:G{0,... ,P} X {0,. ..,T} 

C{x) 

Meshes of different size can be compared after introducing 0(PT) dummy vertices, such 
that the meshes have equal dimensions and each vertex has an aligned vertex on the other 
mesh. For this, it is important to note that any quadrilateral can be subdivided into an 
equivalent mesh of four quadrilaterals with any point on the original quadrilateral as their 
shared corner. 

To extend this further, consider the case where the matching p{p,t) = (’’"(p))''"(f)) is 
defined by two piecewise-linear functions tt and t of |7r| and |t| vertices. This allows comparing 
moving curves under predefined realigments of timestamps as well as positions. For such a 
matching, the surface C{x) = A{x) — B{p{x)) is a quadrilateral mesh of 0{{P + |7r|)(T + |t|)) 
vertices. We illustrate this in Figurej^for f = 4 and |t| = 3 with vertices (t( 0), r(1.5), t( 4)) = 
(0, 2.5,4). In such case 0{T + |r|) timestamps of A are matched with 0{T + |t|) timestamps 


t =( 

1 1.5 2 3 

1 t=0 1 1.5 2 3 4 


\ \ \ \ \ \ \ 

\ \ \ \ N \ \ 

' ' \ \ \ ' ' 

'' \'\ \'\\ '' 


1 1 1 

1 1 1 

1 1 1 

1 1 1 

1 1 1 

1 1 1 

1 1 1 

1 1 1 

1 1 1 

1 1 1 

1 1 1 

1 1 1 

t =( 

1 2 2.5 3 

1 f=6 i 2 2^5 3 4 


B B 

M Figure 3 A piecewise linear matching of timestamps. Left, the matching drawn with the original 
time axes. Right, the matching drawn after warping the time axes. 

The same can be done given a piecewise-linear reparameterization of positions. As a result, 
the Frechet distance under a predefined piecewise-linear reparameterization of timestamps 
and positions can be computed in 0{{P -I- |7r|)(T + |t|)) time. 
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[T] Synchronous Constant Matchings 


In this section, we consider the class of synchronous constant matchings 
where the matching yL{p,t) = (7r(p),t) assumes no realignments of time and 
a constant reparameterizatioiQ tt of positions. Thus, the goal is to find a 
continuous nondecreasing surjection tt on [0, P], such that the Frechet distance 
is minimized. 


Before we present the algorithm to find synchronous constant matchings, 
we refer to an existing algorithm that computes the Frechet distance between static curves. 
This classic algorithm makes use of a data structure called the e-freespace diagram, which is 
the set of parameter pairs for which the represented points are at most e apart. A freespace 
diagram for two static curves is illustrated in Figure]^ 


The matching tt{x) = y between the parameter spaces of static curves A and B can be 
embedded as a bimonotone path {{x,y) \ tt(x) = y} in the freespace. Therefore, the Frechet 
distance is at most e if and only if there exists a bimonotone path from the bottom-left to 
the top-right corner of the e-freespace. Using freespace diagrams, the decision problem for 
static curves can be solved in O(P^) time. 



M Figure 4 Left: Curves A (red) and B (blue) with Frechet distance e. Right: In white, their 
freespace diagram p, = {(a;, y) £ [0, 6] x [0, 5] | |! A(x) — B{y) || < e}. If we draw the freespace diagram 
for any smaller value of e, no bimonotone path through the freespace connects the bottom-left corner 
to the top-right corner. 


We extend this approach to compute matchings between moving curves. For this, consider 
the 3D freespace defined by Equation]^ Since any synchronous constant matching /i(p,t) = 
(7r(p), t) is defined by a bimonotone path tt : [0, P] —>■ [0, P], the matching p is embedded in 
the freespace diagram as the surface p = {{x,y) \ 7r(x) = y} x [0, T], Such a matching yields 
a Frechet distance of at most e if and only if /i C pg. The 3D freespace consists of O(P^T) 
cells Cx,y,t = Pe n ([a:, cc -I- 1] x [y,y+l] x [t, f -|- 1]) for {x, y, t) e N^. 


Pg = {(cr, 2 /,t) G [0,P] X [0,P] X [0,T] | ||A(a;, t) - P(y, t)|| < e} 


(4) 


^ We make slight abuse of notation: we allow y to be drawn from the closure of the space of matchings. As 
a consequence, tt{x) is not an actual function since some x can be matched to multiple (but connected) y 
values. 
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We can simplify the three-dimensional freespace (J>) into a two-dimensional one 
with 


fl {{x,y)G[0,P]x[0,P]\\\A{x,t)-B{y,t)\\<e} (5) 

tG[o.r] 

and find the path tt defining /i in it. To simplify Equation]^ we prove the following lemma. 

► Lemma 1. A cell of the freespace has a convex intersection with any line parallel to 

the xy-plane or the t-axis. 


Proof. Each cell in the freespace is the freespace induced by a pair of quadrilaterals, so 
consider two quadrilaterals A[x, t) and Bly, t). These quadrilaterals are bilinear interpolations 
between four corner points. Hence, A(x,t) — B{yA) is an affine map for each fixed t. 
Likewise, A{x,t) — B{yA) is an affine map for each fixed pair {x,y). Because the preimage 
of a convex norm ball under an affine map is convex, the intersection of a freespace cell with 
lines parallel to the t-axis or the xj/-plane forms a convex set. ◄ 

► Lemma 2. J'P = f] {{x,y) G [0, P] x [0, P] \ \\A{x,t) - B{y,t)\\ < e}. 

Proof. By Lemma the intersection of a cell of with a line parallel to the t-axis is convex. 
Hence, if a point {x,y,t) ^ then either {x,y, [tj) ^ or {x,y, [t]) ^ Therefore the 
internal structure of freespace cells can safely be ignored. ◄ 


We denote the O(P^) cells of the 2D freespace by Cx,y = ([a;, x -|- 1] x [j/, y -|- 1]) n J"™ 
where (x, y) G N^. 

► Lemma 3. Every cell Cx,y in the 2D freespace is convex. 


Proof. By Lemma each cell is the intersection of convex sets, so each Cx^y is convex. ◄ 


Figure illustrates (in white) what the 2D 
freespace B^^ might look like for two moving 
curves. As before, any smaller value of s would dis¬ 
connect the bottom-left from the top-right corner. 
For a given e, for each of the 0{P^) cells, the bound¬ 
ary intervals can be computed in 0{T) time. There¬ 
fore, finding an xy-monotone path takes 0{P^T) 
time, solving the decision problem for the Frechet 
distance under synchronous constant matchings 
in 0{P^T) time. 

3.1 Parametric Search 



M Figure 5 Jy ° for two moving curves. 

To compute the exact Frechet distance efficiently. Darker areas are in fewer layers of 
we use the decision problem in a parametric search 

for the minimum e admitting a matching. When increasing e starting from £ = 0, there are 
three types of critical values of e for which a passage might open in the freespace: 


a) The minimal e with (0,0) G B^^ and {P,P) G B^^. 

b) One of the four boundaries of a cell Cx^y becomes nonempty. 

c) The lower (or left) endpoint on the boundary of cell Cx,y aligns with the upper (or right) 
endpoint on the boundary of Cx+i,y (or Cx,y+j)- 
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Note that all critical values occur when two endpoints align (or an endpoint aligns with a 
gridpoint). For each of the 0{T) layers defining , we use a function in e to represent 
the X- or y-position of such endpoint (or gridpoint). This amounts to a total of 0{P^T) 
functions, and each critical value occurs when two of them intersect. Since the intersections 
between any pair of these functions can be computed in constant time, we can apply a 
parametric search EO]. 

We use Cole’s variant m of the parametric search to find the desired critical value 
in 0{{k + timedec) log k) time. Here k = 0{P^T) is the number of functions to which the 
parametric search is applied and timcdec = 0{P^T) is the running time for the decision 
problem. We obtain the running time of Theorem^ 

► Theorem 4. The synchronous constant Frechet distance between quadrilateral meshes can 
be computed in 0{P^T log{PT)) time. 

► Remark. To compare quadrilateral meshes under a piecewise linear realignment timestamps, 
we can subdivide the quadrilateral meshes as explained in Section Although for simplicity 
we have assumed that the meshes are of equal size, when computing the Frechet distance 
between meshes of different sizes, P x T and Q x T, the 2D freespace has only 0{PQ) 
cells. Thus the decision problem is solved in 0{PQT) time and the exact computation 
takes 0{PQTlog{PQT)) time. 

► Remark. If the inputs of the algorithm are two sequences of 0{T) curves without predefined 
interpolations to obtain quadrilateral meshes, we can still measure their Frechet distance for 
the optimal (but unknown) linear interpolation. Due to the convexity of freespace cells, the 
Frechet distance is the minimum Frechet distance between two curves At and Bt, which can 
be computed in 0{TP^ log P) time by running the original algorithm 0{T) times. 

Synchronous Dynamic Matchings 

Synchronous dynamic matchings align timestamps under the identity match¬ 
ing, but the matching of positions may change continuously over time. Spe¬ 
cifically, the matching is defined as ^(p,t) = Here, ■ 

[0, P] X [0,T] —[0, P] X [0,T] is continuous, and for any t the match¬ 
ing TTt : [0, P] —>■ [0, P] between the two curves at that time is a nondecreasing 
surjection. 

4.1 Freespace Partitions in 2D 

Recall that the freespace diagram is the set pairs of points that are within distance £ of 
each other. 

{x,y) G Pe <S4> \\A{x) - B{y)\\ < e 

If A and B are curves with parameter space [0, P], then their freespace diagram is two- 
dimensional, and the Frechet distance is the minimum value of e for which an cry-monotone 
path (representing p) from (0,0) to (P, P) through the freespace exists. 

We use a variant of the max-flow min-cut duality to determine whether a matching 
through the freespace exists. Before we present the 3D variant for moving curves with 
synchronized timestamps, we illustrate the idea in the fictional 2D freespace of Figure 
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Here, any matching—such as the green path—must 
be an x- and y-monotone path from the bottom left 
to the top right corner and this matching must avoid 
all obstacles (i.e. all points not in J'e). Therefore each 
such matching divides the obstacles in two sets: those 
above, and those below the matching. 

Suppose we now draw a directed edge from an 
obstacle a to an obstacle b if and only if any matching 
that goes over a must necessarily go over b. The key 
observation is that a matching exists unless such edges 
can form a path from the lower-right boundary to the 
upper-left boundary of the freespace. In the example, 
a few trivial edges are drawn in black and gray. If all 
obstacles were slightly larger, an edge could connect a 
blue obstacle with a red obstacle, connecting the two 
boundaries by the edges drawn in black. 



\ZK 


M Figure 6 A matching (green) in the 
2D freespace (white). 


4.2 Freespace Partitions in 3D 


In contrast to the 2D freespace where the matching is a path, matchings of the form /i(p, t) = 
{nt{p),t) form surfaces in the 3D freespace (see Equation]^. Such a surface again divides 
the obstacles in the freespace in two sets and can be punctured by a path connecting two 
boundaries. We formalize this concept for the 3D freespace and give an algorithm for deciding 
the existence of a matching. 


{x,y,t) G Te if and only if \\A{x,t) — B{y,t)\\ < £ 


( 6 ) 


For x,y,t & N, the cell Cx,y,t of the 3D freespace is the set J>n([a;, a;-|-l] x [y, y+1] x [t, t+l])- 
The property of Lemma [l] holds for all such cells. 

We divide the set of points not in Bg into a set O of ^ ^ 

so-called obstacles, such that each individual obstacle is 
a connected point set. Let u be the open set of points 
representing the left and top boundary of J-g. Symmetrically, 
let d represent the bottom and right boundary, see Figure]^ 

Denote by O' C O the obstacles between the boundaries. 

0 = {u,d}U0' with U O' = ([0,P]2 X [0,r])\ J-,; 
u = {{x,y,t) I (x < 0 A y > 0) V (x < P A y > P)}; 
d = {{x, y,t) I (x > 0 A y < 0) V (x > P A y < P)}. 

Given a matching /i, let P C O be the set of obstacles 
below it, then u ^ D and d G D. Here, we use axes (x, y, t) 
and say that a point is below some other point if it has a smaller y-coordinate. Because each 
obstacle is a connected set and y cannot intersect obstacles, a single obstacle cannot lie on 
both sides of the same matching. Because all matchings have u ^ D and d G D, a matching 
exists if and only if ^{d G D ^ u £ D). 

We compute a relation > of elementary dependencies between obstacles, such that its 
transitive closure © has d © u if and only iidGD^uGD. Let a > 5 if and only if a U 6 
is connected (a touches b) or there exists some point {xa,ya,ta) G a and {xb,ybAb) G b 
with Xa < Xb, ya > yb and ta = tb- We prove in Lemmasandthat this choice of > satisfies 
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the required properties and in Theorem that we can use the transitive closure © of > to 
solve the decision problem of the Frechet distance. 

► Lemma 5. If a © &, then a G D ^ b G D. 

Proof. Assume that a>b, then either a touches b and no matching can separate them, or 
there exists some {xa,ya,t) G a and {xb,yb,t) G b with Xa < Xb, ya > yb- If there were 
some matching y with a G D, then (xa,ybnt) G y for some > ya- Similarly, if & ^ D, 
then (xb,y'^,t) G y for some y'^ < yb- We can further deduce from Xa < Xb and monotonicity 
of y that we can pick y'^ such that ya < y^i ^ y'^ < yb- However, this contradicts ya ^ yb, 
so such a matching does not exist. Hence, a G D ^ b G D whenever a> b and therefore 
whenever a © b- ◄ 

► Lemma 6. If d G D ^ u G D, then d©u- 

Proof. Suppose d G D ^ u G D but not d © u- Then no matching exists, and no path 
from d to u exists in the directed graph G = {0,t>). Pick as D the set of obstacles reachable 
from d in G, then D does not contain u. Pick the tightest matching y such that D lies below 
it, we define y in terms of matchings C x {t} in the plane at each timestamp t- 

{x, y, t) G TTt if and only if (x > x Ay' < y) ^ ^m{x', y', t) A m{x, y, t) where 

m{x,y,t) if and only if {{x' ,y' ,t) \ x' < x Ay' > y} C\\^D = % 

Because u ^ D, this defines a monotone path ttj from (0, 0) to (P, P) at each timestamp t. 
Suppose that tt^ properly intersects some o GO, such that some point of {xo,yo,t) G o lies 
below TTt- It follows from the definition of > and —•m{xo,yo,t) that d>o for some d G D- 
However, such obstacle o cannot exist because D satisfies >. As a result, no path intersects 
any obstacle and we can connect the paths TTt to obtain a continuous matching y without 
intersecting any obstacles. So y does not intersect obstacles in O \ P, contradicting d G D ^ 
uG D- ◄ 

► Theorem 7. The Frechet distance is greater than e if and only if d © u for e. 

Proof. We have for every matching that u ^ D and d G D- Therefore it follows from 
Lemma that no matching exists if d © u for e. In that case, the Frechet distance is greater 
than e. Conversely, if ^(d © u) there is a set D satisfying © with u ^ D and d G D- In that 
case, a matching exists by Lemma and the Frechet distance is less than e. ◄ 

We choose the set of obstacles O' such that [jO' = ([0, P]^ x [0, P]) \ and the relation t> is 
easily computable. Note that due to Lemmaeach connected component contains a corner 
of a cell, so any cell in the freespace contains constantly many (up to eight) components 
of U O' - Moreover, we can index each obstacle in O' by a grid point {x, y, t) G N^. 

Let Ox^y,t C ([0, P]^ X [0, T]) n ([x — 1, a; + 1] x [y — 1, y + 1] x [t — 1, t + 1]) \ be the 
maximal connected subset of the cells adjacent to (a;,y,t), such that contains {x,y,t)- 
Now, the obstacle o^^y^t is not well-defined if {x,y,t) G P^, in which case we define o^^y^t to 
be an empty (dummy) obstacle. We have O' = ^ remark that obstacles 

are not necessarily disjoint. 

Each of the 0{P^T) obstacles is now defined by a constant number of vertices. We 
therefore assume that for each pair of obstacles (a, 6) G O^, we can decide in constant time 
whether a>b; even though this decision procedure depends on the chosen distance metric. 
For each obstacle a = Ox,y,t, there are 0{P‘^) obstacles b = Ox',y',t' for which a>b, namely 
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because t — 2<t'<t + 2iia\>b. Furthermore, u and d contribute to 0{P^T) elements of 
the relation. Therefore we can compute the relation > in 0{P'^T) time. 

Testing whether d ® u is equivalent to testing whether a path from d to u exists in the 
directed graph (0,[>), which can be decided using a depth first search. So we can solve the 
decision problem for the Frechet distance in 0{P^T + | > |) = 0{P'^T) time. However, the 
relation > may yield many unnecessary edges. In Section [4.4| we show that a smaller set E of 
size 0{P^T) with the same transitive closure © is computable in 0{P^T log P) time, so the 
decision algorithm takes only 0{P^T log P) time. 


4.3 Parametric Search 


To give an idea of what the 3D freespace looks like, we have 
drawn the obstacles of the eight cells of the freespace between 
two quadrilateral meshes of size PxT = 2x2inF igure 
Cells of the 3D freespace lie within cubes, having six faces and 
twelve edges. We call such edges x-, y- or f-edges, depending 
on the axis to which they are parallel. 

We are looking for the minimum value of e for which a 
matching exists. When increasing the value of e, the relation > 
becomes sparser since obstacles shrink. Critical values of e 
occur when > changes. Due to Lemma ^ all critical values 
involve an edge or an xt-face or j/t-face of a cell, but never the 
internal volume, so the following critical values cover all cases. 



■ Figure 8 [0,2]^ \ 


a) The minimal £ such that (0,0, t) G Te and (P, P, t) G for all t. 

b) An edge of Cx,y,t becomes nonempty. 

c) Endpoints of y-edges of Cx,y,t and Cx+i,y,t align in y-coordinate, or endpoints of x-edges 
of Cx,y,t and Cx,y-j,t align in x-coordinate. 

d) Endpoints of a t-edge of Cx,y,t and a t-edge of Cx+i,y-j,t align in t-coordinate. 

e) An obstacle in Cx,y,t stops overlapping with an obstacle in Cx+i,y,t or Cx,y-j,t when 
projected along the x- or y-axis. 

The endpoints involved in the critical values of type a), b), c) and d) can be captured 
in O(P^r) functions. We apply a parametric search for the minimum critical value £abcd 
of type a), b), c) or d) for which a matching exists. This takes 0{{P^T + timcdec) log(Pr)) 
time. 

We illustrate the need for critical values of type e) in Eigurej^ 
here obstacle a overlaps with both obstacles b and c while the 
overlap in edges does not contribute to >. It is unclear how critical 
values of type e) can be incorporated in the parametric search 
directly. Instead, we enumerate and sort the 0{P^T) critical values 
of type e) in 0{P^Tlog{PT)) time. Using 0(log(PP)) calls to the 
decision algorithm, we apply a binary search to find the minimum 
critical value £e of type e) for which a matching exists. Finding 
the critical value £e then takes 0((P^T + fe'medec) log(Pr)) time. 

The synchronous dynamic Frechet distance is then the minimum 
of £abcd and £e. This results in the following running time. 

► Theorem 8. The synchronous dynamic Frechet distance can be computed in 0{{P^T + 
ijmcdec) log(PT)) time. 

Before stating the final running time, we present a faster algorithm for the decision algorithm. 
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4.4 A Faster Decision Algorithm 


To speed up the decision procedure we distinguish the cases for which two obstacles may be 
related by >, these cases correspond to the five types of critical values of Section |4.3| Critical 
values of type a) and b) depend on obstacles in single cells, so there are at most 0{P^T) 
elements of > arising from type a) and b). Critical values of type c) and e) arise from pairs 
of obstacles in cells in the same row or column, so there are at most 0{P^T) of them. In 
fact, we can enumerate the edges of type a), b), c), and e) of t> in 0{P^T) time. On the 
other hand, edges of type d) arise between two cells with the same value of t, so there can 
be 0{P‘^T) of them. 

We compute a smaller directed graph {V,E) with \E\ = 0{P^T) that has a path from d 
to u if and only if d © m. Let V = O = {u, d} U O' be the vertices as before (we will include 
dummy obstacles for grid points in that lie in the freespace) and transfer the edges in [> 
except those of type d) to the smaller set of edges E. We must still induce edges of type d) 
in E, but instead of adding 0{P'^T) edges, we use only O(P^T) edges. The edges of type d) 
can actually be captured in the transitive closure of E using only 0{P) edges per obstacle 
in E. 

Using an edge from Ox,y,t to Ox+i,y,t and to Ox,y-i,t, we construct a path from Ox,y,t to 
any obstacle Ox+i,y-j,t- The sole purpose of the dummy obstacles is to construct these paths 
effectively. For obstacles whose gridpoints have the same t-coordinates, it then takes a total 
of 0{P^T) edges to include the obstacles overlapping in t-coordinate related by type d), this 
is valid because {x,y,t) G Ox,y,t for non-dummy obstacles. 

Denote by the edges of type d) of the form (a, b) = {ox,y,ta i Ox+i,y-j,th) where tb = ta+k, 
then the set Eq of O(P^T) edges is the one we just constructed. Now it remains to induce 
paths with ta ^ h, that still overlap in t-coordinates, i.e. the sets E'^ 2 t and E^. 

Denote by t~(a) and the minimum and maximum t-coordinate over points in an 

obstacle a. For each obstacle, both the t~{a) and the t~^{a) coordinates are an endpoint of a 
t-edge in a cell defining the obstacle due to Lemma ^ and thus computable in constant time. 

Our savings arise from the fact that Eq induces a path 
from Ox+i^y—j ^t+k to Ox-\-i\y—j' if Ox^y^t ^ ^x+i^y—i^t+k 
and Ox,y,t > Ox+i',y-j',t+k with i < i' and j < f, so 
we do not need an additional edge to induce a path 
to the latter obstacle. To avoid degenerate cases, we 
start by exhaustively enumerating edges of E^ {k G 
{—2, —1,1,2}) for which z < 1 or j < 1 in 0{P^T) time 
so we need only consider edges with z > 2 A j > 2. 

For these remaining cases, we have a > 5 if and only 
if t'^{a) > t~{b) Atb = ta + k, and t~{a) < t^{b) A 
tb = ta — k for positive k. From this we can derive the 
edges of E'^. Although for each a, there may be O(P^) 
obstacles b such that a>b with tb = ta + k, the Pareto 
frontier of those obstacles b contains only 0{P) obstacles, 
see the grid of fictional values t~{b) in Figure 10 In 
Section |4.5| we show how to find these Pareto frontiers 
in 0(P log P) time per obstacle a, using only 0{P'^T) 
preprocessing time for the complete freespace. 

As a result, we can compute all 0{P^T) edges of in O(P^riogP) time. By Theorem]^ 
the decision problem for the synchronous dynamic Frechet distance is solvable in 0{P^T log P) 
time. 


y 



H Figure 10 Two edges (green) 
cover (red) all four obstacles h (green) 
within the query rectangle (blue) 
with values t~ {b) < t^(a) = 2. 
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► Theorem 9. The decision problem for the synchronous dynamic Frechet distance is solvable 
in 0{P^T log P) time. 

Proof. The edges E of types other than d) are enumerated in 0{P^T) time, and using 
constantly many Pareto frontier queries for each obstacle, 0{P^T) edges of type d) in E are 
computed in 0{P^TlogP) time. Given the set E of edges, deciding whether a path between 
two vertices exists takes 0(|i?|) = 0{P^T) time. The transitive closure of E equals ®, so a 
path from d to u exists in E if and only if there was such a path in [>. Since we compute E 
in 0{P^T log P) time, the decision problem is solved in 0{P^TlogP) time. ◄ 

The following immediately follows from Theorems and 

► Corollary 10. The synchronous dynamic Frechet distance can be computed in time 
O{P^Tlog Plog{PT)). 


4.5 Pareto Frontier Queries 

Suppose we are given a matrix M with m columns and n rows, and we want to be efficiently 
query submatrices for the Pareto frontier of numbers that are at most a given threshold value, t. 
A query specifies the threshold t, and two coordinates {xo,yo) and {xi,yi) of the query 
rectangle R = {ccq, ..., Xi} x {yo ,..., yi}. Let Ct(R) = {(a;, y) G i? | M[x, y] < t} denote the 
coordinates of the query rectangle that must be dominated by the Pareto frontier Ft (i?) C 
Ct{R). That is, if {x,y) G Ct{R), then some {x' ,y') G Ft{R) with x' < x Ay' <y exists. 

We preprocess each of the 0{n) rows of M in 0{m) time by storing their cells as the 
leaves of an augmented binary tree, whose internal nodes store the minimum value over 
its subtrees. Then queries for the index of the leftmost element with value at most t in a 
range {xq, ... ,a;i} of that row can be answered in 0{logm) time. We can compute Ft by 
including for each row, the element with minimum a;-coordinate and value at most t in the 
query range. So using 0(n) queries, we compute a set Ft{R) of size 0(n) in 0{nlogm) time, 
using 0(nm) preprocessing time for the matrix M. 

In this case, Ft{R) is not actually the Pareto frontier since some of its elements might 
be dominated by other elements. With a slight modification, we can make Ft{R) the 
actual Pareto frontier (of minimum size). If the query with range {xq, ... ,xi} returns x* 
with xo < X* < xi for some row, then the query range for subsequent rows can be restricted 
to {xo ,..., a;* — 1}, so no unnecessary values are generated. 


5 


Hardness 


We extend the synchronous constant and synchronous dynamic 
classes of matchings (of Sections and to asynchronous 
ones. For this, we allow realignments of timestamps, giving 
rise to the asynchronous constant and asynchronous dynamic 
classes of matchings. The asynchronous constant class ranges 
over matchings of the form = (7r(p),T(<)) where the tt 

and T are matchings of positions and timestamps. The asynchronous dynamic class of 
matchings has the form = (7rt(p), r(t)) for which the positional matching irt changes 

over time. We first prove that the asynchronous constant Frechet distance is in NP. 



► Theorem 11. Computing the Frechet distance is in NP for the asynchronous constant class 
of matchings. 
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Proof. Given any matching /i(p, t) = (7r(j>), r(t)) with a Frechet distance of e, we can derive— 
due to Lemmajl]— a piecewise-linear matching t* in 0{T) time, such that a matching t) = 
(tt* (p) , T* {t )) with Frechet distance at most e exists. We can realign the quadrilateral meshes A 
and B under r* to obtain meshes A* and B* of polynomial size. Now the polynomial-time 
decision algorithm for synchronous constant matchings (see full paper) is applicable to A* 
and B*. ◄ 


Due to critical values of type e), it is unclear whether each asynchronous dynamic matching 
admits a piecewise-linear matching r* of polynomial size, which would mean that the 
asynchronous dynamic Frechet distance is also in NP. 

We show that computing the Frechet distance is NP-hard for both classes by a reduction 
from 3-SAT. The idea behind the construction is illustrated in the two height maps of 
Figure 11 These represent quadrilateral meshes embedded in and correspond to a single 


clause of a 3-CNF formula of four variables. 



U Figure 11 Two quadrilateral meshes A and B embedded in (indicated by color and isolines). 
Their Frechet distance is 2 isolines if the clause {X 2 V - 1 X 3 V 1 X 4 ) is satisfiable and 3 isolines 
otherwise. 




M Figure 12 The freespace of (A, B) at times (0,0) (left) and (T/2,T/4) (right). 


We distinguish valleys (dark), peaks (white on A, yellow on B) and ridges (denoted W, Fi 
and Ti). An important observation is that in order to obtain a low Frechet distance of e < 3, 
the n-th valley of A must be matched with the n-th valley of B. Moreover, each ridge Xi 
must be matched with Fi or Ti and each peak of A must be matched to a peak of B. Note 
that even for asynchronous dynamic matchings, if Xi is matched to Fi, it cannot be matched 
to Ti and vice-versa because the (red) valley separating Fi and Ti has distance 3 from Xi. 

The aforementioned properties are reflected more clearly in the 2D freespace between 
the curves at aligned timestamps t and r(t). In Figure 12, we give two 2D slices (with 


{IaAb) = (0,0) and (TI2,TjT), respectively) of the 4D-freespace diagram with £ = 2 for the 
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shown quadrilateral meshes. In this diagram with e = 2, only 2^ monotone paths exist (up 
to directed homotopy) whereas for e = 3 there would be 2^ monotone paths (one for each 
assignment of variables). For £ = 2, the peak of X 2 cannot be matched to F 2 at t = T/4 
of B, corresponding to an assignment of X 2 = true. 

Consider a 3-CNF formula with n variables and m clauses, then A and B consist of m 
clauses along the t-axis and n variables {Xi ... and Fi,Ti... Fn,Tn) along the p-axis. 
The fc-th clause of A is matched to the fc-th clause of B due to the elevation pattern on the 
far left (p = 0). This means that the peaks of A are matched with peaks of the same clause 
on B and all these peaks have the same timestamp because r(t) is constant (independent 
of p). 

For each clause, there are three rows (timestamps) of B with peaks on the ridges. On each 
such timestamp, exactly one ridge (depending on the disjuncts of the clause) does not have a 
peak. Specifically, if a clause has Xi or -^Xi as its fc-th disjunct, then the fc-th row of that 
clause has no peak on ridge Fi or Ti, respectively. We use these properties in Theorem [T^ 
where we prove that it is NP-hard to approximate the Frechet distance within a factor 1.5. 

► Lemma 12. The Frechet distance between two such moving curves is at least 3 if the 
corresponding 3-CNF formula is unsatisfiable. 

Proof. Consider a matching yielding a Frechet distance smaller than 3 given an unsatisfiable 
formula, then the peaks of A (of the fc-th clause) are matched with peaks of B (of a single row 
of the fc-th clause). Assign the value true to variable Xi if ridge Xi is matched with Ti and 
false if it is matched with Fj. Then for every clause (Vi V F,- V 14) with Vj G {Xi,^Xi}, there 
is a peak at 7r(Ai), TT{Xj) or 7r(Afe) for that clause. Such a matching cannot exist because 
then the 3-CNF formula would be satisfiable, so the Frechet distance is at least 3. ◄ 

► Lemma 13. The Frechet distance between two such moving curves is at most 2 if the 
corresponding 3-CNF formula is satisfiable. 

Proof. Consider a satisfying assignment to the 3-CNF formula. Match Xi with the center 
of Fi or Ti, if Xi is false or true, respectively. For every clause, the timestamp with peaks 
of A can be matched with a row of peaks on B. As was already hinted at by Figure El the 
remaining parts of the curves can be matched with £ = 2. Therefore this yields a Frechet 
distance of at most 2. ◄ 

► Theorem 14. It is NP-hard to approximate the asynchronous constant or asynchronous 
dynamic Frechet distance for moving curves in within a factor 1.5. 

Proof. By Lemmas and [T^ the asynchronous constant or asynchronous dynamic Fre¬ 
chet distance between two quadrilateral meshes embedded in is at least 3 or at most 2, 
depending on whether a 3-CNF formula is satisfiable. ◄ 


5.1 Orientation-Preserving Homeomorphisms 


Previous results HI [ID] have shown that computing the Frechet distance 
between surfaces under orientation-preserving homeomorphisms is NP-hard 
for surfaces embedded in K^. We will refer to this variant as the orientation¬ 
preserving Frechet distance. The prior results hold for triangular meshes, 
which are a degenerate case of quadrilateral meshes. Although the orientation¬ 
preserving Frechet distance is upper-semicomputable, it is unknown whether 
it is computable at all. We consider the case where the corners of the meshes are matched 
with each other. 






K. Buchin, T. Ophelders and B. Speckmann 


15 


The prior NP-hardness constructions for the orientation-preserving Frechet distance seem 
unnecessarily complex. Therefore we extend our results for the asynchronous dynamic Frechet 
distance of Section to obtain a new hardness construction for the orientation-preserving Fre¬ 
chet distance for surfaces embedded in Note that we cannot directly apply the previous 
construction because there was only a single matching of timestamps for asynchronous 
dynamic matchings. Thus, we do not preserve the property that for a clause, the timestamp 
with peaks of A maps to a single timestamp of 13. 

We prevent this by means of a second dimension in which we enforce that a row of peaks 
of A maps to a single row of peaks of B. In addition to the embedding in of Figure 11 


which defines one coordinate for every point on a clause of a quadrilateral mesh, we define a 
second coordinate using Figure 13 yielding an embedding in 



H Figure 13 The second coordinate of a clause gadget of A (top) and B (bottom) embedded in R^, 
the first coordinate is given in Figure 1111 


Now, under the maximum norm a row of peaks of A can only be matched to 

a single row of peaks on B for e < 3. Conversely, we can still match the meshes of two 
satisfiable formulas with e = 2. Hence, Theorem [^follows. This result extends to triangular 
meshes since all quadrilaterals lie in the plane, and can thus be represented by a pair of 
triangles. For norms other than the maximum norm, the problem is still NP-hard, but our 
bound on the approximation factor is smaller than 1.5. 


► Theorem 15. Unless P=NP, no polynomial time algorithm can approximate the orientation¬ 
preserving Frechet distance between two quadrilateral meshes embedded in under the 
maximum norm within a factor 1.5. 


6 


Conclusion 


Based on the Frechet distance, we presented several similarity measures between moving 
curves, together with efficient algorithms for computing some measures, while proving 
NP-hardness for computing others. 

Although many algorithmic solutions to variants of the Frechet distance between static 
curves also apply to the synchronous constant Frechet distance, extending the synchronous 
dynamic class is more complex. For example, computing the synchronous dynamic Frechet 
distance between moving closed curves (with a cylinder instead of square as parameter space) 
remains an open problem. 

Finally, the performance of the algorithm for the synchronous dynamic Frechet distance 
seems undesirably slow and we have not been able to prove tight lower bounds on this 
running time yet. 
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